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Abstract 

We study Skyrmion dynamics in a (2+l)-dimensional Skyrme model. The system 
contains a dimensionless parameter a, with a = corresponding to the 0(3) sigma- 
model. If two Skyrmions coUide head-on, then they can either coalesce or scatter — 
this depends on a and on the incident speed v, and is affected by transfer of energy to 
and from the internal vibrational modes of the Skyrmions. We classify these internal 
modes and compute their spectrum, for a range of values of a. In particular, we find 
that there is a fractal-like structure of scattering windows, analogous to those seen 
for kink-antikink scattering in dimensions. 
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1 Introduction 



In this paper, we investigate aspects of the dynamics of the planar Skyrme system, in 
which the field configuration is a map from to S^. In particular, we examine the 
discrete vibrational modes of static A^-Skyrniion solutions. Such modes correspond 
to relatively long-lived vibrational states. They are of importance for semiclassical 
quantization (see, for example, [19, 2] for the Skyrme case, [20] for the planar-Skyrme 
case); but our interest lies in their effect on the classical dynamics of Skyrmions. 

A long-standing prototype is the system in 1 dimensions. Here the soliton 
is a kink (or antikink), which possesses a single internal oscillatory mode. If a kink 
and an antikink approach each other with low relative speed v < vq, then they form a 
long-lived 'breather' or bion state [9]. (This is not an exact breather, and eventually 
decays.) If, on the other hand, the kink and the antikink approach at high speed 
V > vi, then they bounce off each other, and each escapes to infinity (the collision 
is inelastic, with some radiation being emitted). For intermediate impact speed 
vo < V < vi, there is a fractal-like structure of 'reflection windows', with trapping and 
reflection alternating [4, 3, 1]. This can be understood in terms of a resonant energy 
exchange between the translational motion of the kink and its internal oscillation [4] . 
Also worth mentioning is that the internal mode can in some sense be 'extrapolated' 
to a non-infinitesimal dynamical process, namely a kink-antikink-kink collision [14]. 

In this paper, we shall study analogous features for the (2+l)-dimensional Skyrme 
system. There are several significant differences between this case and that of kinks. 
The first is that two kinks (as opposed to a kink and an antikink) always repel each 
other, and there is no static 2-kink solution; whereas for Skyrmions (with a suitable 
choice of potential) static 2-Skyrmion solutions do exist. Consequently, it makes sense 
to study Skyrmion-Skyrmion collisions: for head-on collisions, one might expect that 
there will be a critical speed vq such that 

• for impact speed v < vq, the Skyrmions coalesce to form a 2-Skyrmion; 

• foi V > Vq, the Skyrmions scatter and each escapes to infinity; 

and the internal modes should play a role in this process. 

As we shall see, however, the picture is not quite as simple as this. The system 
contains a dimensionless parameter a, and the various dynamical features depend 
crucially on a. We study numerically how the spectrum of vibrational modes, and 
the scattering behaviour, vary with a. For small a, internal modes are absent, and 
the picture is indeed as suggested above: coalescence for v < vq, and scattering 
for V > Vq. But for larger a, a rich spectrum of internal modes appears, and the 
scattering behaviour also becomes more complex. In particular, there is a range of 
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a within which one sees a fractal-hke structure of 'scattering windows', separated by 
regions of coalescence. 



2 The planar Skyrme system 

The Skyrme system in R^+^ is defined as follows. Let x'* = (x^ ^x^ ^x^^ = {t,x,y) 
denote the standard space-time coordinates; indices are raised and lowered using the 

standard Minkowski metric (with signature H ). The two spatial coordinates are 

denoted . The Skyrme field is a unit vector field (f) = {(pi, (f>2, ^s), with (j) ■ (j) = 1. 
Its space-time and spatial derivatives are denoted d^cp and djcp respectively. The 
Lagrangian density is 

c = \{d^$) ■ {d^$) - 11^,,,^^"' - \aV{4>^) , (1) 

where Q,^^ is the triple scalar product Vt^y = (j) ■ [d^cp) x {dy(j)), y is a function of 03, 

and a and 7 are constants. 

The boundary condition at spatial infinity \s (f) ^ (0, 0, 1) as r — > 00, where 
= + y^. A configuration satisfying the boundary condition has an integer 

winding number, which we denote N, and which is given by 

iV = ^ y O12 dPx . (2) 

The static energy £^ of a configuration is 

E = \j [{dj$) ■ {dj$) + 7(^^12)' + aVicP^)] (fx . (3) 

The quantity (7/a)^/"^ has units of length, and so we shall henceforth fix the length 
scale by taking 7 = 0. So we have a system depending on the parameter a, as well 
as on the function V{4>^). The static energy of a configuration with winding number 
N satisfies the Bogomol'nyi bound [7, 6] 



£; > 47r|iV| ( 1 + I y" v^F(^#) . (4) 



Clearly for finite energy we need V{1) = Q. In the asymptotic region r ^ 1, the 
two components cpi and ^2 (which are the analogues of the three pion fields in the 
full Skyrme model) satisfy a Klein-Gordon equation where the 'pion mass' m is given 
by rrf = —aV'{l)/2. So the frequency u of radiation is bounded below by m. Two 
choices of V{(f)) for which the corresponding systems have been investigated in some 
detail [17, 18, 11, 21, 16] are V{(p) = 1 - (/> and V{4)) = 1-^2. j.g£gj. ^-^ese as 
Old Baby Skyrme (OBS) and New Baby Skyrme (NBS), respectively [21]. We shall 
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restrict our attention to these two systems, concentrating especially on the NBS case; 
note that previous work on semiclassical quantization [20] dealt with the OBS case. 

We remark in passing on the a ^ oo limit, which amounts to deleting the (5/^0)^ 
term in the Lagrangian [15]. This system may admit stable solitons. For example, if 
V{(j)) = 1 — (j), then there is a solution which has compact support {ie. = 1 outside 
a disc of finite radius), and which docs not saturate the Bogomol'nyi bound [5]. With 
V{(f)) = (1 — 0)^, on the other hand, one has a smooth static rotationally-symmctric 
solution with = 1 — 2exp(— r^/2A''), which does saturate the Bogomol'nyi bound. 
See [15] for a more general discussion. 

The simplest Skyrmion solutions are rotationally-symmetric, or more accurately 
0(2)-symmetric. Letting r and 9 represent the usual polar coordinates on R^, we 
say that a configuration (t>{r,9) is 0(2)-symmetric if 

• <l>z = 4>2.{r) with ^3(00) = 1 = 1^63(0)!; and 

• + i(j)2 = F{r) cxp(iA^0) with F{r) real-valued. 

Here N is an integer. Note that wc necessarily have F{oo) = = F[0). An equivalent 
definition is to say that has the 'hedgehog' form 

${r, e) = (sin(/) cos(Ar^), sin(/) sin(Ar^), cos(/)) , (5) 

where the profile function f = f{r) is smooth and real- valued with /(oo) = and 
/(O) = Kit for some positive integer K. There exist symmetric solutions with K > \, 
but they are unstable [11]; so we shall restrict our attention here to the case K = 1. 
In this case, N is the same as the winding number (2) [17]. For the NBS system, 
there is considerable numerical evidence that, for each N, there is a smooth minimal- 
energy iV-Skyrmion solution, and this solution is 0(2)-symmetric. For a = 1, the 
normalized energy = E/{AttN) of this A^-Skyrmion for 1 < iV < 4, obtained by 
numerical minimization, is as follows: Ei = 2.15, E2 = 1.91, -E3 = 1.85, £'4 = 1.83. 
(Note that the Bogomol'nyi bound (4) is E^ > 1.7854.) Since E^ is a decreasing 
function of A^, we expect Skyrmions to coalesce; in particular, a low-speed collision 
of two l-Skyrmions will result in a single 2-Skyrmion. 

3 Vibrational modes of the A^-Skyrmion 

In this section, we study the spectrum of vibrations about 0(2)-symmetric Skyrmion 
solutions. The set of all perturbations about a symmetric configuration is a vector 
space which is acted on by 0(2), and which therefore decomposes into disjoint sub- 
spaces corresponding to representations of 0(2). We begin, therefore, with a brief 
summary of the irreducible representations of 0(2) in this context. 
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The group 0(2) is generated by the rotations i— > 6+c, which make up S0(2); and 
the reflections 9 h- > —0. It has two 1-dimensional irreducible representations (+) and 
(— ); and for each positive integer p, a 2-dimensional irreducible representation (p), 
containing two degenerate modes (p, +) and {p, —). The corresponding perturbations 
5(j) of (j) have the following form. 

(+): 5$= (Acos(iV0),^sin(iV0),5), with Asin/ + Bcos/ = 0; 
(-): 5$= {-Asu^{Ne),Acos{Ne),{))■ 
(p,+): 5^i=Bcos[{N -p)e)] + Ccos[{N + p)e)], 

S(f>2 = B sin[(Ar - p)e)] + C sin[(iV + p)e)] , 

S(f)3 = A cos{pe), with {B + C) sin / + ^ cos / = 0. 

(p,-): 5(t)i = -Bsm[{N -p)9)]-Csm[{N+p)9)], 
6(j)2 = B cos[{N - p)e)] + C cos[{N + p)9)], 
S(f)^ = -A sm{pe) , with {B + C) sin / + ^ cos / = 0. 

Clearly (+) consists of perturbations which maintain the hedgehog form (5): the 
proflle of the Skyrmion changes, but it remains 0(2)-symmetric. And (— ) consists of 
perturbations which allow F{r) to become complex- valued but preserve its modulus 
(in other words, F acquires an r-dependent phase). So (— ) includes the zero-mode 
01+102 ^ exp (ix) (01 +1^2)7 where x is a constant phase angle. The other zero-modes 
are the translations in space; these form a doublet which belongs to the representation 
(1). Apart from shape modes in (+), some other particularly significant positive 
modes (cf. [19, 2] for the (3+l)-dimensional Skyrme case) are 'dipole breather' modes 
belonging to (1), which correspond to the Skyrmion oscillating from one side to the 
other; and the 'splitting mode' of the 2-Skyrmion into two single Skyrmions, which 
belongs to (2). 

The spectrum of oscillations around a static solution will consist of the zero- 
modes mentioned above (angular frequency uj = 0); a continuum of radiation modes 
(a; > m); a finite number of negative modes (uj imaginary) if the solution is unstable; 
and a finite, possibly zero, number of discrete positive modes {0 < uj < m). 

We have used two, very different, numerical methods to compute the vibration 
spectrum. The first involves deriving the relevant Sturm-Liouville equation (in r) for 
each of the classes (+), (— ), (1), (2), . . .; this is then solved using a Chebyshev spectral 
method. This method can be used for unstable as well as stable solutions. The second 
method involves a full (2-1-1 )-dimensional simulation of the field equations, and only 
works for vibrations about stable solutions. The procedure is to construct a static 
Skyrmion (by relaxation), put in a perturbation, solve the time evolution for a long 
time interval, extract a time-series by sampling the field at some point in space, and 
finally Fourier-analyse this data. See, for example, ref [2], where this method was 
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used for the (3+l)-dimensional Skyrme system. The results of the two methods are 
consistent. 

Let us now look at the internal modes of the 1-Skyrmion in both the OBS and 

NBS systems. The results arc summarized in Figure 1, where the modes can be seen 
'peeling off' the bottom of the contimium band as a increases. For the OBS system 
(where = a/2), there arc no internal modes if a < 1. The first mode to appear 
belongs to (1), followed by modes in (2), (— ) and (+) in that order. For the NBS 
system (where = a), the first mode (again in (1)) appears at a « 0.27, followed 
by a mode in (+) (at a ^ 0.30). The latter crosses over the former at a 0.31; 
but at a 1.2, the modes cross again, and the (1) mode has the lower frequency for 
a > 1.2. The number of modes grows quite rapidly with a. 




In the OBS system, the 2-Skyrmion is stable; for all o; > 0, it has a discrete 
positive mode in (2), corresponding to the breakup of the 2-Skyrmion into two single 
Skyrmions ('splitting mode'). As a increases, more positive modes appear. The 3- 
Skyrmion in the OBS system is unstable [17]; there is always a negative mode, which 
belongs to (2). From now on, we consider the NBS system only. 

For the NBS system with a = 1, the positive discrete modes of the 1-Skyrmion 
and the 2-Skyrmion are as follows. 

A'^ = 1: a shape mode (in (+)) with uj = 0.806; and a dipolc breather doublet in (in 
(1)) with UJ = 0.834. See the right-hand diagram in Figure 1. 

A'^ = 2: a shape mode with u = 0.519; a dipole breather doublet with u = 0.613; 
a splitting doublet with lo = 0.398; and two further doublets (one in (2) with 
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u = 0.909 and one in (3) with u = 0.912). See the right-hand diagram in Figure 
3. 



4 Skyrmion-Skyrmion scattering 

In the NBS system, sohtons attract each other, and the lowest-energy N = 2 configu- 
ration corresponds to two superimposed sohtons. This means that if two sohtons are 
released from rest at finite separation, they slowly move towards each other to form 
a large lump. When simulating this numerically, we see that after the scattering, the 
sohtons emerge at 90°, but their mutual attraction slows them down and forces them 
to move back towards each other. They then merge again and scatter in their initial 
direction, but stop again and keep oscillating back and forth. As the interaction is 
non-elastic, the lumps emit waves that carry away some energy. As a result, the 
amplitude of oscillation slowly decreases, and the configuration converges toward the 
two-soliton static configuration. 

If the solitons are sent toward each other with a large-enough speed, they mo- 
mentarily merge to form a large lump, and then two lumps emerge at 90° and move 
away from each other. If the initial speed is too small, then these solitons do not 
escape from their mutual attraction, and they eventually form a single lump after 
several oscillations. 

This behaviour suggests that there is a critical speed below which the solitons 
merge, and above which they escape. In Figure 2, we present the dependence of the 
critical speed v (defined below) on the coupling constant a. In the limit a — > 0, 
the system is the 0(3) sigma model, and a head-on collision always results in 90° 
scattering [13, 12]. In other words, we expect that — > as a ^ 0, and this is 
consistent with the numerical results. For small a, the Skyrmions scatter very easily; 
but as a increases towards 0.27, they have an increasing tendency to coalesce. In a 
collision between Skyrmions, any positive discrete (localized) mode of the Skyrmion 
will inevitably be excited (cf. [3, 8]); this removes translational kinetic energy from 
each Skyrmion, and transfers it to an internal excitation which eventually decays into 
radiation. So when the shape mode appears (at a = 0.27), it very efficiently removes 
translational energy from the Skyrmions, and leaves them in a bound state. There 
is a partial effect even below a = 0.27, which can perhaps be understood in terms of 
a quasimode (cf. [8]): an 'internal' mode which is embedded within the continuum 
spectrum, and which becomes the shape mode when a increases beyond 0.27. 

On the other hand, for a larger than about 0.7, the critical speed decreases. 
Between these two extremes, there is a region where the critical speed is close to 
the speed of light (1 in our units); but as scattering two solitons numerically at 
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large speed is quite difficult, we were not able to determine the critical speed when 
0.27 < a < 0.7, or determine for which value of a the 90° scattering never takes 
place. 



0-8 I I I I I I I I I I I I I I I I I I I I I I I I I 




0.1 — 

I I I I I I I I I I I I I I I I I I I I I I I I I I 
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Figure 2: Critical velocity as a function of a for the NBS model. 



If one looks closely at Figure 2, one notices that the curve folds back around the 
value a = 0.263. In that region, the solitons merge when the initial speed is either 
below the lowest value or above the largest one. The numerical simulations from 
which Figure 2 was obtained involve the head-on collision of two Skyrmions in their 
most attractive mutual orientation (the 'attractive channel'), and with each having 
equal initial speed. The critical speed v plotted in Figure 2 is defined as the speed 
below which the solitons always merge, except for the top part of the folding where 
it is defined as the speed above which the solitons always merge. 

When a is slightly smaller than 0.263, and when the initial speed is slightly above 
the lower critical value, we observed numerically that the scattering behaviour of the 
solitons can have a complex pattern for some values of a. Instead of simply escaping 
or merging, the solitons swing back and forth as if they were slowly merging, but 
eventually they escape at 90° or in the direction from which they came. 

To describe this process, we can define the 'scattering number' as the number 
of times the two solitons scatter at 90°, and define it as zero when the solitons 
coalesce. So when the initial speed is very large, the scattering number is 1 (except 
for the narrow range of a around 0.26, where high-speed solitons coalesce). When 
the scattering number is even, the solitons scatter in the direction they came from; 
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and when it is odd, they scatter at 90°. For a range of values of a, the scattering 
number is very sensitive to the initial velocity, and exhibits a fractal-like structure. 

The structure of the data is described in Table 1, where we present the scattering 
number for various ranges of initial speed, sampled at regular intervals. Each digit 
in the large strings corresponds to the scattering number for a given value of v. The 
final string (which is split over two lines) should be read as a single merged string. 
The first digit in each string corresponds to the smallest value of v ('First in the 
table), while the last digit corresponds to the largest value of v ('Last f ' in the table). 
The increment size in v is given in the third column. 

One clearly sees that for the values of a given in Table 1, there is no sharp 
transition between scattering and the merging of solitons, but that instead there 
are several windows of initial speed where the solitons scatter several times before 
escaping at 90° or in their original directions; these windows are separated by regions 
where the solitons do not escape. As the scattering number increases, the width of 
the windows become smaller, but the structure at any scale seems to be similar. 
One can also see in Table 1 that often, at the edge of a given window, there is a 
narrower window for a larger scattering number; this is not always the case, but we 
have observed this many times. 

We should also point out that the scattering patterns given in Table 1 form a small 
subset of the range of scattering velocities that we have investigated numerically. 
We have observed a fractal structure in the scattering data for the range of speed 
contained between the lower and upper branch of the critical speed on Figure 2 when 
a was just below 0.26. When a = 0.24, the windows were not as rich in structure 
as for a = 0.26, exhibiting large gaps of coalescence between the regions of multiple 
scattering. 

5 Scattering and Vibration Modes 

Our observations are analogous to those described by Anninos et al [1] in their study 
of kink-antikink scattering in the (p^ model in one dimension, where the scattering 
also exhibits windows of scattering separated by regions where a long-lived bound 
state is formed. 

It is interesting to note that in Figure 2, the value a = 0.263 where the critical 
speed curls over is just below the value a = 0.265 where the breather mode (1) 
emerges from the continuum band. Moreover, the curve for the shape mode (+) 
emerges at a = 0.295, and the curves for the two modes cross at a = 0.310. See 
Figure 3. 

We also notice from Figure 3 that a splitting mode in (2) (although not the lowest 
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one) peels off from the continuum band at a = 0.62, and this is approximately the 
value (cf Figure 2) above which the solitons can scatter, suggesting that this splitting 
mode is implicated in the large-a part of Figure 2. 

To understand the relation between the vibration modes and the scattering of 
solitons, we must first realise that during the scattering process the solitons are ex- 
cited in various vibration modes. When they are sufficiently far apart, each lump 
can be considered as a single Skyrmion; while when they overlap, they form a con- 
figuration close to the N = 2 rotationally-symmetric solution. Between these two 
configurations, they form intermediate states that extrapolate between the two ex- 
treme configurations, but that are difficult to describe. One thing that can be inferred 
from the expressions for (50 is that the mode for the 1-Skyrmion becomes the 
mode {2p)2 for the 2-Skyrmion, as the two Skyrmions merge. For example, the dipole 
breather mode for each incident 1-Skyrmion becomes the splitting mode (2)2 for 
the 2-Skyrmion, as the Skyrmions coalesce. More precisely, the lowest state for 
which is the translation mode of the 1-Skyrmion becomes the lowest-state of the (2)2 
mode. The first excited state of the (l)i mode, shown in Figures 1 and 3, transforms 
into the first excited state for the (2)2 mode, which (as shown in Figure 3) crosses 
the mass threshold around 0.62. 

Moreover, due to the nonlinear nature of the model, the different vibration modes 
exchange some energy during the scattering. In particular, this means that some 
kinetic energy is transfered from the translation mode into the various vibration 
modes, and the reverse is also true when the two Skyrmions try to escape from their 
mutual attraction. 
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During the scattering, the vibration modes that are excited vary with time and 
depend on the overlap between the Skyrmions; but when the initial speed is close to 
the critical value, we observed that the vibration period is less than the characteristic 
time of scattering, implying that the system has plenty of time to oscillate. 

We should also point out that all the scatterings that wc have simulated nu- 
merically involved head-on collisions in the attractive channel [18], which means in 
particular that (j) was invariant under the rotation 9 9 + tt. This implies that the 
modes (p) of the 2-Skyrmion, with p odd, were never excited. 

We thus see that when a < 0.265, there is no internal mode to excite the 1- 
Skyrmion, and the oscillation of the Skyrmions induced by their scattering is radi- 
ated away. When a > 0.265, the is excited during the scattering and it can 
absorb progressively larger amounts of energy. As the Skyrmions merge, this mode 
transforms into the first exited (2)2 mode which is well above the mass threshold. 
The energy is thus radiated away and the scattering docs not take place. When 
a > 0.65, the (2)2 is a bound state, hence the energy is not radiated away and the 
90° scattering becomes possible again. 

To try to predict the critical velocity shown in Figure 2, wc now consider a 
simple model for the two-soliton scattering process, namely the dynamics of two 
equal masses, connected by a spring, on a cone. The cone corresponds to the moduli 
space for the two Skyrmions in their symmetric configuration [10], while the two 
masses connected by the spring models a single internal mode of oscillation of a 
Skyrmion. Moreover, to model the attraction between the two Skyrmions, the two 
masses are made to evolve in a central potential, describing the binding energy of 
two Skyrmions. Each mass is equal to half the mass of a Skyrmion, and the masses 
are positioned so that the line joining them crosses the origin, ie the tip of the cone. 
In this picture, the Skyrmion scatters with its mirror image; and the move over the 
tip of the cone corresponds to a 90° scattering once the cone is unfolded onto the 
plane. 

The Lagrangian for this system is 



The constants A and A are coefficients, depending on a, which we will choose to 
model the two-soliton scattering as closely as possible. The parameter Xq describes 




(6) 



where 



(7) 
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the size of the soliton; we will take Xq = 1. Notice that the depth of V2 is A/2; 
so given that there are two degrees of freedom, A corresponds to the depth of the 
two-Skyrmion bound-state potential. The equations of motion are 

AAsinh(AXi) 



= -k{Xr 
M^X2 = k{Xi-X2) 



2cosh2(AXi) ' 
.4Asinh(AX2) 
2cosh2(AX2) ■ 



(8) 



The value of k is chosen so that the frequency of oscillation between Xi and X2 

is the frequency of the shape mode (+), namely k/M^ = a;<+>j. The equation 
a;<-|.>-^ = 0.30 + 0.63(1 — exp(— 1.58 a)) gives a good approximation. Note that 
our simple model does not distinguish between the different vibration modes; as 
the frequencies of the (+) (1) modes are very similar, we have taken the former for 
convenience. 

The parameters A and A in the potential V2 are determined as follows. For A 
we take the depth of the potential, namely A = 1 — E2/{2Ei), where E2 is the 
energy of the two-soliton bound state and Ei is the energy of a single soliton. In 
Figure 4 we show the a-dependence of the binding energy E2/{2Ei). The curve is 
well-approximated by the relation E2/{2Ei) = 1 + 0.12(exp(-3.88 a) - 1) + 0.01 a. 
For A, we impose the condition that the frequency of the small-amplitude oscillations 



_ 0.95 




Figure 4: a-dependence of the binding energy of two Skyrmions. 



for V2 is equal to the frequency of the lowest splitting mode for the two-soliton 
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bound state a;<2>2 ■ That frequency is well approximated by the expression a;<2>2 = 
0.2 + 0.2(1 - exp(-3.875 a)) - 0.015a. We therefore have A = a;<2>2(2MtotM)^/^ 
where Mtot = 2Mr. 

To be more realistic, we should have one oscillator for each vibration mode. One 
quick way to approximate this is to have no decoupled oscillators which, to first 
approximation, have the same frequency. This will correspond to taking 2no points 
of mass Mx/uq linked in pairs by springs with elastic constant k/uQ. To preserve 
the total binding energy of the system we must also divide A by no. This is actually 

1/2 

equivalent to solving equation (8), after multiplying A by Uq . 

To simulate a scattering, we set up the two masses so that they are separated by 
their equilibrium distance Xq, and so that their centre of mass is located at x = 10. 
We then send both of them with the same speed towards the origin. The motion 
of the masses in the potential well stretches the string and results in some transfer 
of translation energy into the oscillator. The initial speed therefore has to be large 
enough for the two masses to go over the tip of the cone and escape towards x = — oo. 
If the speed is too small, the two masses oscillate around the tip of the cone. 

Using these parameters, we have determined the critical velocity for no = 1 and 
no = 2, as shown in Figure 5. It produces the correct shape of curve, ie one that 
looks like the inverse of a Morse potential, but the actual critical velocities are too 
small. This is explained by the fact our simple model does not take into account the 

0.3 
0.25 
0.2 

^ 0.15 
0.1 
0.05 


0.5 1 1.5 2 2.5 

a 

Figure 5: Critical speed for the 2-point model; one oscillator (bottom curve) and two 
oscillators (top curve) 
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radiation of energy. Otherwise, the shape of the curve is more or less explained by 
the shape of the binding energy of two Skyrmions, as shown in Figure 4. Where the 
well is deepest, around a = 1, more energy is transferred into the oscillator, and the 
critical speed is high. 

The drop in the critical velocity for 0.27 < a < 0.3 when no = 2 is caused by a 
phase resonance between the oscillation of the system in the potential well and the 
oscillation of the rigid oscillator. Similar phenomena were also observed when taking 
a different value for Xq. 

Figure 1 suggests that our toy model predicts the critical speed reasonably well 
for large values of a, if we take no somewhere between 1 and 2. For small values of a, 
the critical speed is too small by roughly a factor of 2. The predicted critical speed 
is also far too small in the range 0.3 < a < 0.7, but the position of the maximum is 
surprisingly in the correct region of a. 

The main success of our simple model is in explaining how the existence of a 
critical speed for 90° scattering comes from the fact that some kinetic energy is 
transferred into the vibration modes of the system. It also shows that the dependence 
of the critical speed on a is related to the depth of the potential well between the two 
Skyrmions. To explain the other features of the curve shown on Figure 2, one must 
analyse how the vibration modes for 1-Skyrmions and 2-Skyrmions transform into 
one another, and consider when these modes are above or below the mass threshold. 

When a < 0.265, our simple model does not really apply, as the 1-Skyrmion 
does not have any genuine vibration mode, although the 2-Skyrmion does. When 
0.265 < a < 0.62, the the mode transforms into the (2)2 which radiates its 
energy away and makes the 90° scattering difficult or impossible. 

When a > 0.62, none of the major modes excited during the scattering radiates 
energy away. So some of the energy stored in these vibration modes can be converted 
back into the translation mode, and the 90° scattering can happen at a relatively 
small speed. 

To predict the critical speed more accurately, one would have to take into account 
all the modes that are excited in the process, as well as how they are coupled together, 
and coupled to the deformation of the system during the scattering. This goes well 
beyond what we can expect from such a simple model. One could consider using 
a genuine geodesic approximation for the Skyrme model, but as we do not have an 
analytic expression for the general 2-Skyrmions configuration with a separation, this 
is difficult to do. 

The fractal structure of the scattering data just below a = 0.265 is more difficult 
to explain. When the vibration modes have a frequency just above the mass thresh- 
old, the solution for the continuum spectrum exists for all frequencies, but for some of 
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them, the eigenfunction has a pronounced maximum at the position of the Skyrmion. 
These solutions can be thought of as oscillation modes that decay relatively slowly 
in the linear limit. These quasi-modes can thus absorb some energy and make the 
system oscillate. The fractal structure of the scattering data probably comes from a 
delicate phase resonance between these modes and the scattering oscillations. 

6 Conclusion 

We have shown that the planar Skyrme model possesses some interesting vibration 
and scattering properties. In the limit of small-amplitude oscillations, the Skymions 
can vibrate at a finite number of frequencies that depend on the coupling constant 
a. 

When two Skyrmion collide head-on, they can coalesce, scatter at 90° or, in some 
very special circumstances, scatter back-to-back. For small values of a, the Skyrmions 
can scatter at 90° if their collision speed is large enough. When 0.265 < a < 0.62 the 
90° scattering speed seems impossible, within the limit of our numerical integration, 
while for a > 0.62 the 90° scattering becomes possible again. When a is just below 
0.265, the scattering data has a fractal structure. 

We explained the possibility of 90° scattering by studying the vibration modes 
of the 1-Skyrmion and 2-Skyrmions solutions. When the deformation mode of the 
1-Skyrmion transform into the excited splitting mode of the 2-Skyrmion, and when 
this mode is above the mass threshold, a large amount of energy is radiated and the 
90° scattering is not possible. We also described a simple dynamical model which 
shows that the value of the critical speed is related to the binding energy of the 
two-Skyrmion solution. 
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a 


First V 


V increment 


Last V 


Scattering Numbers 


0.24 


0.4705 


0.0001 
002222000000( 


0.477 
3002220050 


00022000003220000200002002002030200000111111 


a 

0.24 


First V 
0.4721 


V increment 
0.00001 


Last V 
0.4726 
2222222 


Scattering Numbers 
22222222220000000000000000500000555500000000 


a 

0.24 


First V 
0.472428 
0555555 


V increment 
0.000001 
555000000000( 


Last V 
0.4725 
)000000000 


Scattering Numbers 
00000000000000000000700000000000055555555555 


a 

0.25 


First V 
0.49 

22222 


V increment 
0.0001 
200000000002: 


Last V 
0.497 
2222000000 


Scattering Numbers 
00222250000022250000222000020002200200202002 


a 

0.25 
2222^ 


First V 
0.493241 
>000000000 


V increment 
0.000001 
000000005000( 


Last V 
0.49332 
3000005000 


Scattering Numbers 
00005000000000005505555500000040055043333000 


a 

0.26 


First V 
0.4905 
00002 


V increment 
0.0005 
022203002000( 


Last V 
0.53 
3022222000 


Scattering Numbers 
00002220000022230002200002000220230202022200 


q; 

0.26 


First V 
0.5156 


V increment 
0.000005 

2222222^ 


Last V 
0.5159 
>222222222 


Scattering Numbers 
22222222222220000000003000300000003330000000 


a 

0.26 


First V 
0.51578 
0300 


V increment 
0.000001 
000000000333^ 


Last V 
0.515872 
1000000000 


Scattering Numbers 

00000000333333040000000000000000000000000000 

00005333333333333334040 



Table 1: Scattering numbers for a range of a and v. 
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